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ABSTRACT 

We have used a Lagrangian, fully implicit, one-dimensional 
hydrodynamic computer code tc investigate the evolution of a gas 
cloud impacting the surface of a 20 km, 1 M q neutron star. This 

cas is initially at rest with respect to the surface of the neutron 
star, extends to 185 km above the surface, and is optically thick. 
The infall results in a burst which lasts about 0.1 seconds and 

c 

reaches a peak luminosity and effective temDerature of 2.4 x 10“ L„ 

6 0 
and 9 x 10 K; respectively. The burst was followed by a phase of 

oscillations with a period 0.2 seconds. 

INTRODUCTION 

One of the more interesting scenarios for the cause of the X- 

ray Bursts is the magnetospheric gate theory of Lamb, et al_ J In 
this scenario we have a rapidly rotating neutron star with a strong 
magnetic field in a binary system with a companion that is losing 
mass toward the neutron star. This material is impeded by the mag- 
netic field in such a fashion that accretion is not continuous but 
occurs in episodic events. To our knowledge there have been no 
previous numerical simulations of this phenomena. ~ 

We use the same hydrodynamic code described in the last paper . 
Because it is a Lagrangian, stellar interior, computer code which 
assumes radiation transport by diffusion; we cannot place a hole in 
the mesh and so assume that the blob of gas extenas from the surface 
of the neutron star to a radius of ^ £5 km. This blob is oDtically 

thick, has a mass of 10"^ ** , and is initially at rest with respect 


to the surface of the neutron star. The mass of this blot was cho- 
sen so as to simulate a much smaller amount of material falling 
only onto the polar caps. The gas initially has a temperature of 

% 10 y K and a surface luminosity of 0.1 L . We begin the calcula- 
tion by allowing the material to fall onto the surface. 

RESULTS 

Once the infall has begun, an accretion shock forms at the 
boundary of the core. After ahout lOOu sec of evolution, the in- 
fall velocities have reached free fall and the temperature behind 

g 

the sno.ck has reached 10 K. It takes 3.5 milliseconds for the 
outer edge of the shock to reach the surface of the blob which has 
now fallen inward to 80 km. The entire envelooe bounces imparting 

-14 

enough energy to the two outermost zones (2 x 10 M ) for them to 

reach escape velocity. The light curve for the burst is shown in 
Figure l . The peak luminosity and effective temperature produced 

* 5 6 

by the infalling material is 2.4 x 10 L Q and 9 x 10 K, respective- 
ly. Although, this effective temperature is too low for any of the 
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observed bursts, we could raise this temoerature by reducing the 
inner radius of the neutron star and reducing the amount of inf ail- 
ing material. 

The terr.Derature in the shell source at the core-envelope inten- 
ts 

face (CEI ) has remained at 10 K and enough triole-a reactions are 

16 

occurring to keep the rate of nuclear energy generation at 10 erg 

gm sec . In addition, the deeper layers of the envelope are 
convective, which is mixing fresh unburned nuclei into the shell 
source which also is maintaining the large rate of nuclear burning. 
However, the outer layers are expanding too rapidly for the energy 
produced in the shell source tc reach them and their luminosity 
crocs rapidly. For a short period the s^-decays in the surface 
layers are able to halt the decreasing luminosity but the rapid ex- 
pansion again overcomes this energy source. The entire burst lasts 
about 0.3 sec. 

The decline in luminosity continues until the escaDing shells 
become optically thin. At this time the energy from the interior 
can now escape and the luminosity begins to increase. Now that we 
can see into the deeper layers we find that the accreted material 
remaining on the neutron star is pulsating with a period of '■« 0.2 
sec (Figure 2). The average luminosity of this object is --10 4 * L 

and the effective temperature is 0.3 kev. It has a radius of 150km 
(average) . 

These pulsations are a direct result of the infall and are net 
caused by a partial ionization mechanism such as in the Cepheias or 
ZZ Ceti variables. We also find that these pulsations must be 
occurring in high overtones since we find a number of nodes in the 
eigenfunction and the largest amplitudes are in the surface zones. 

We closely followed these pulsations for 4 seconds of star 
time. However, the computer time was becoming prohibitive and we 
ended this part of the calculation. We removed the escaping zones 
and followed the resulting evolution for another thousand seconds. 
No further changes occurred and we ended the calculation. It is 
estimated that at least one day of further evolution would be nec- 
essary to completely burn the hydrogen in the envelope. 

We also studied a purely cooling sequenc® in which we followed 
the contraction of the extended envelope assu... ng only neutrino 
losses - no nuclear energy production. The envelope contracts from 


about 203 km tc 13C km in 1C J jec with the luminosity climbing to 2 


x 1C 5 L and the effective temoerature to about 1 kev. It takes 

G 

another 4C minutes to contract to 20 km. During the cooling chase 
the surface luminosity reaches a peak luminosity of 2 x 10 3 L Q and 


a peak effective temperature of 1.9 kev. 


CONCLUSIONS 


We find that our Burst was too cool and too rapid to resemble 



the noma! bursts although the tine scale is in agreement with 

some of tne rapic bursts that have oeen observed for some objects'. 
It is certainly possible tnat more realistic initial conditions, 
such as assuming the material was already infall i ng , could im- 
prove the correspondence between cur calculations and the observa- 
tions. 

We also found that the infall produced a ringing of the enve- 
lope which showed short period pulsations. Such pulsations have 

been observed in at least one burst event 3 . 
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